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Coupling dark energy to dark matter provides one of the simplest way to effectively modify gra- 
vity at large scales without strong constraints from local (i.e. solar system) observations. Models 
of coupled dark energy have been studied several times in the past and are already significantly 
constrained by cosmic microwave background experiments. In this paper we estimate the constraints 
that future large-scale observations will be able to put on the coupling and in general on all the 
parameters of the model. We combine cosmic microwave background, tomographic weak lensing, 
redshift distortions and power spectrum probes. We show that next-generation observations can 
^ ' improve the current constraint on the coupling to dark matter by two orders of magnitude; this 

O . constraint is complementary to the current solar-system bounds on a coupling to baryons. 
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I. INTRODUCTION 



• Several observational campaigns will be launched in the next fewyears to advance our knowledge of the matter 
■ distribution from redshifts of order z < 0.5 (the depth of the SDSS [l| and 2dF surveys) to < 2 — 3 (the range 
of satellite missions like Euclid [3| and of ground-based surveys like LSST [|, DES 0, etc). These surveys will map 
hundreds of millions of galaxy redshifts and billions of galaxy images distributed across almost the full accessible sky. 



The information contained in these data will occupy the cosmologists for years to come and it is very useful already at 
O ' this stage to foresee what these data can tell us about the currently investigated cosmological models. Several works 
^ ' have been devoted to forecasting future constraints with future wide surveys and a small selection of the recent ones 
^ ■ can be found in Ref. jft-ilOj. 

' The main science drivers for most of these large-scale enterprises is the quest for dark energy (DE) . Even if we get 
very accurate knowledge of the cosmic background expansion and therefore of the dark energy equation of state, this 
will not be enough to distinguish among competing models. In fact, in the models in which dark energy clusters or 
^ , gravity is modified (see e.g. [ll|). the linear growth of clustering is in general not completely fixed by background 
' observables. Mapping the linear clustering has been identified as the main probe of modification of gravity or, more 
[ in general, of non-standard dark energy models. 

The main probes of matter clustering are the weak lensing and the galaxy power spectrum P{k) and indeed most 
large scale surveys aim at mapping one or both of these quantities. For instance Euclid, probably the most ambitious 
project to date, sets as goal the mapping of half the sky both in imaging and in broad-band spectroscopy, with an 
average redshift of order unity and a maximal useful redshift around 2 or 3. Combining weak lensing and P(fc) gives 
the opportunity to break several degeneracies among the cosmological parameters and to directly measure the growth 
function. 

Both WL and P{k) probe the universe at relatively small redshifts. A much longer lever arm is obtained by 
1 combining them with the cosmic microwave background, in particular employing the specifics of the Planck mission. 
?-j I In this way we can tighten sensibly the constraints on all the cosmological parameters. 

In this paper we focus on one particular class of dark energy models that includes typical features of modified 
gravity models, the coupled dark energy model ( (T2| - [l^ . see e.g. [ll| for a review). The idea behind this model is 
very simple: when dark energy is promoted to a scalar field 0, as in the vast majority of DE models and also as in 
the prototypical accelerated expansion, inflation, one naturally wonders whether this field could be mediating a force 
among matter particles. This is indeed the same phenomenology as in the Brans-Dicke model, except that now we 
require this field also to drive acceleration and therefore to be dominant and to carry a potential. In this paper we 
model the potential as an inverse power-law V ~ as in the Ratra-Peebles scenario [ISj . This one-parametric 
family is simple but gives us enough freedom to explore the dependence of the results on the potential slope and 
therefore on the equation of state. The phenomenology of this interaction can be immediately grasped by writing 
down the energy-momentum conservation equations; 

TL).:,-(^T^^)^-^^ (1) 



T^r^).:. = -/^TW'/';. (2) 
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where the subscripts m, (p denote matter and dark e nergy, r espectively and where /3 is the couphng constant. While 
/3 could also be a coupling function /3((/)) (see e.g. |16l - [l8| ). in this paper we restrict ourselves for simplicity to a 
constant value. A coupled particle of mass m will respond to this interaction by following in the Newtonian regime 
the scalar-gravitational potential 

Cm 

r 

where depends on the field potential. For a dark energy field one expects the interaction scale l/m^ to be of 
astrophysical size. The parameter quantifies therefore the deviation from Newtonian dynamics. Any indication 
that /3 ^ would of course show a fundamental modification of nature's law. The main goal of this paper is to forecast 
the future constraints on /? and to see how all the other cosmological parameters depend on it. 

Since gravity now is supplemented by an additional scalar mediator, the resulting scalar-tensor force deviates from 
Einstein's (and from Newton's) one and this gives rise to a number of observable differences with respect to the 
uncoupled models, both at the background and at the perturbation level. The expansion is generally modified when /3 
is non-zero. In particular, it has been shown in [l^ that for most potentials DE is not negligible in the past, contrary 
to the cosmological constant case and to models that do not drastically deviate from it. Adopting the Ratra-Peebles 
potential one finds that the evolution of the scalar field during the matter domination is controlled by the coupling 
/3 and the density fraction $7^ is proportional to As a consequence, the expansion rate H, given in terms of the 
effective equation of state Wcg by 

|^ = -^(l + ^^cff) (4) 

is modified. We define here Wcft = Pt/Pti where pt and pt are the total pressure and energy density contributions 
respectively. The value of Wef / decreases from 1 /3 in radiation dominated era to 

Was ~ ^/3' (5) 

during matter dominated epoch, where the dominant contribution to pt is given by the scalar field. More recently, 
dark energy dominates and wde tends to —1: the higher the coupling, the fastest is the transition between matter 
and dark energy dominated era. This behavior is illustrated in Fig.([T]) Accordingly, properties as the age of the 
Universe and the distance to last scattering (and therefore the position of the acoustic peaks) do also depend on (3 
and this allows to put stronger constraints on it. In fact the best constraints on /3 so far come from the CMB data of 
WMAP, < 0.15 [2O, '2ll|. We stress that here we consider massless neutrinos but it is important to notice that if 
neutrinos are massive higher values of the coupling (/3 < 0.2) can be allowed [22]. At the perturbation level, the extra 
pull induced by the scalar force adds to gravity and speeds up the perturbation growth, resulting in further possible 
constraints. 

Any modification of gravity has to pass the strong bounds from local experiments, both in laboratory and in the 
solar system or other astrophysical sources (see e.g. (23l - [25| ). This can be achieved in several ways by the so-called 
chameleon (26l . [27| or Vainshtein mechanisms [l^ but here we adopt a simpler way out: we assume that baryons 
are uncoupled. In this case, all local experiments are bypassed and cosmology becomes the only way to observe the 
coupling. 

II. COUPLED DARK ENERGY 

We consider coupled dark energy cosmologies as described by the lagrangian: 

1 
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C = ^-d^cf>d^cf> - C/(0) - m{<f>)^^ + C^inm , (6) 



in which the mass of matter fields tp coupled to the DE is a function of the scalar field (j>. In the following we will 
consider the case in which DE is only coupled to cold dark matter (CDM, hereafter denoted with a subscript c). The 
choice m(0) specifies the coupling and as a consequence the source term Q(0)^ via the expression: 

dlnmjct)) 

<9(0)m = Pc 5^0- (7) 

Due to the constraint of conservation of the total energy-momentum tensor, if no other species is involved in the 
coupling, = 
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Figure 1: Evolution of for various values of /3. 



The conservation equations for the energy densities of each coupled species are: 

p'^ = -37^^0(1 + - (3(0)0 , (8) 
p'c = -3Hpc + Q(4,)Q 

(where primes denote differentiation with respect to conformal time and H — aH is the conformal Hubble function) 
plus the standard conservation equation for baryons. Here we have treated each component as a fluid with T'^ [a)ii = 
{Pa -\~ Pa)u^u'^ +Pq,(5^, where = (—a, 0,0,0) is the fluid 4-velocity and Wa = Pa/ Pa is the equation of state. The 
class of models considered here corresponds to the choice: 

m(0) = rnoe~^^ , (9) 

with the coupling term equal to 

<3(0)o = -^Pc0' ■ (10) 
Equivalently. the scalar field evolves according to the Klein-Gordon equation: 

0" + 2H0' + a2^ = a2^p, . (11) 
dcf) 

Throughout this paper we choose an inverse power law potential defined as: 

V = ^or" (12) 

with a and Vb constants. When the coupling is absent, this potential leads to a well-known transient tracking solution 
D^, JS,, ^1 in which 

2 

WDE ~ — T (13) 

a + I 

which sets up just before the final dark energy domination; after this regime, asymptotically wde ^ ^1- When 
/3 ^ but small, this tracking regime is also approximately true but squeezed between the modified matter era with 
Weff ~ 2/3^/3 and the final phase with wde — ^ —1 and therefore barely visible. Current supernovae la constraints 
impose an upper limit a < 0.2. 

Various choices of couplings have been investigated in literature [1243 l30l - [36| . Analysis of the models and con- 
straint on these couplings have been obtained in several ways, inclu ding spherical collapse ([s^, Ull and references 
therein), time renormalizazion group A^-body simulations (itI. |40|. |41|| and effects on supernovae, CMB and cross- 
correlation of CMB and LSS [7i. ,a. il9n2a li^-liil. Refs. 0,[i| are particularly similar in spirit to our work. However 
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there are several differences. First, the couphng is different since we use the couphng induced by a scalar-tensor 
model. Second, we include the galaxy clustering probe (baryon acoustic oscillations or BAO, redshift distortions and 
full P{k) shape). Third, we include the potential slope a as an additional parameter. Fourth, we use the updated 
Euclid Definition phase [3] specification (there is some change with respect to the ones in the Assessment phase defined 
in the yellow book in particular, the total area has been reduced from 20,000 to 15,000 square degrees). 

III. METHODOLOGY 

In this section we present the methodology we used to derive predictions on cosmological parameter constraining 
power from a combination of upcoming and future datasets, namely galaxy power spectrum sliced in several redshift 
bins, CMB angular power spectra (temperature and polarization) and Weak Lensing power spectrum. The statistical 
method chosen to do forecasts is the Fisher Matrix analysis. In order to derive the contour confidence regions in the 
parameter space both the likelihood of the data and the distribution of parameter values is assumed to be Gaussian. 

Suppose to have an observable Xobs = [^i, ^jv] which is a function of the parameter set © = [^i, On]). Then 
the probability of estimating unknown parameters based on known outcomes is 

£(Xobs/e) oc exp (Xobs - X(0))^C-\Xobs - X(0)) (14) 

where X(0) is the theoretical prediction for Xobs as a function of 9 and C is the covariance matrix of the observed 
components. By assuming priors on each parameter, one is in principle able to sample this function at reasonable 
accuracy level. 

In the Fisher Matrix formalism the observed outcome is the mean values of the observables X^t assumed as the null 
hypothesis. This method allows a quick way to estimate errors on cosmological parameters, given errors in observable 
quantities. The Fisher matrix is defined as the Hessian of the log-likelihood function £, 

such that if the parameters are assumed to be gaussianly distributed their likelihood can be written as 

£(0)cxexp-i^0,i^.,0,. (16) 

ij 

By the Cramer-Rao inequality, a model parameter 6i cannot be measured to a precision better than (Fu)^^^'^ when 
all other parameters are fixed, or a precision {{F~^)iiY^'^ when all other parameters are marginalized over. In practice, 
the Fisher matrix is a good approximation to the uncertainties as long as the likelihood can be approximated by a 
Gaussian, which is generally the case near the peak of the likelihood and therefore in cases when the parameters are 
measured with small errors. Conversely, if the errors are large, then the likelihood is typically non-Gaussian, and the 
constraint region is no longer elliptical but characteristically banana-shaped. In this case, the Fisher matrix typically 
underestimates the true parameter errors and degeneracies, and one should employ the full likelihood calculation 
approach to error estimation. 

In this work the observables X correspond to the galaxy power spectrum measured by a future experiment like Eu- 
clid, the imminent Planck CMB angular power spectra (TT, EE, BB and TE components) and the weak lensing power 
spectrum as detected by a future Euclid-like mission. The cosmological parameter set 6 comprises the dimensionless 
Hubble parameter h, the dark matter and baryon density parameters flc and fib-, respectively, the spectral index of 
primordial perturbations n^, the coupling constant /3 and the slope of the dark energy potential a. We assume a flat 
geometry. In the Fisher formalism, all the information about the experiment is embedded in the covariance matrix: 
hence we need to define the suitable covariance for each observable/experiment. 

As a fiducial model we assume the values of parameters shown in TablH Other fiducial parameters needed for the 
various probes are specified later on. Since the main cosmological effects depend on (for small (3) we use rather 
than /3, as parameter. This confines the parameter volume to > 0. Choosing the fiducial at = 0, in the Fisher 
approximation the cut > does not alter the confidence regions. Also, we choose a = 0.2 as fiducial, but this 
should not have an important impact since the tracking regime lasts very shortly; we tested also the case a = 0.1. 
The other fiducial values are taken from WMAP7 although since a ^ our fiducial model is not exactly ACDM. 

All the derivatives of the Fisher matrix are performed numerically by evaluating the spectra at two values 9i{l ± e) 
where di is a cosmological parameter (for parameters whose fiducial is zero we use di zt e). Here we chose e = 0.03 
but we also tested that for e = 0.06 the flnal constraints change by at most 20%. 
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Parameter 


Value 




0.1116 


Ho 


70.3 km s"^ Mpc"^ 




0.02227 


Us 


0.966 







a 


0.2 


A 


2.42 X 10"® 



Table I: Cosmological fiducial model: parameters are consistent with the WMAP 7 year results for a ACDM cosmology [4^ . |47| 
plus a coupling. 



Planck 


Channel [GHz] 


FWHM 


GT[^lK/K] 


/sky = 0.85 


70 


14' 


4.7 




100 


10' 


2.5 




143 


7.1' 


2.2 



Table II: Planck experimental specifications 481. FWHM is the Full Width at Half Maximum of the beam, assuming a Gaussian 
profile and expressed in arc-minutes. The polarization sensitivity is ap — 

IV. RESULTS FROM CMB DATA 

The coupling has two main effects on the CMB: 1) it moves the position of the acous tic p eaks to larger P& due to 
the increase in the last scattering surface distance (sometimes called projection effect, [1J| and references therein); 
2) it reduces the ratio of baryons to dark matter at decoupling with respect to its present value, since coupled dark 
matter dilute faster than in an uncoupled model. Both effects are clearly visible in Fig. ^ for some values of /3. 

We use the Fisher Matrix method described in section IIIII to predict confidence contours for the cosmological 
parameter set O = {(3^ ,a,Q,c,h,Q,i,,ns,\ogA\, where A is the spectrum normalization, around the chosen fiducial 
model. The fiducial value of A is chosen to be ^ 2.42 x 10"*^, as in WMAP7 best fit for ACDM OS^. For our 
CMB analysis, we use Planck satellite fis'] specifications and include three frequency channels, each characterized by 
the experimental specifications illustrated in TablTTl For each channel we specify the FWHM (Full Width at Half 
Maximum) of the beam when assuming a Gaussian profile and the temperature sensitivity ut- The polarization 
sensitivity is given by crp = ^/2aT- To each theoretical C/ spectrum, we add a noise spectrum given by: 

Ni^{eafe^^[l{l + l)/lb') (17) 

where l^, is given by lb = y/WhPl/e and 9 is the FWHM. 

We obtained the CMB spectra by modifying the CAMB code [49*1 including the coupling; the output have been 
compared to an independent code [2Q] that is build on CMBFAST and the agreement was better than 1%. The initial 
conditions needed to obtain the desired present values of the cosmological parameters must be found by trial and 
error, through an iterative routine. 

The main results are in Fig. ([3]), obtained marginalizing over all parameters except those in the axes. Notice the 
sharp degeneracy 0^ vs. /i, in agreement with [20|, and with Vic ans Ug. This can be understood in the context of 
coupled dark energy. With respect to a cosmological constant model, coupled dark matter dilutes more rapidly so 
that p„i was higher in the past, leading to a faster expansion and a consequent smaller size of the sound horizon at the 
last scattering surface, not fully compensated by the faster expansion after decoupling; this moves the acoustic peak 
towards larger £s, an effect which is degenerated with an increase of h (Fig. ([3]), central panel). On the other hand, a 
stronger coupling decreases the peak amplitudes because pb/ Pm at decoupling gets smaller: this effect is compensated 
by lower values of Vic or by higher values of the spectral index Ug. In Table (jIVp we report the fully marginalized la 
errors for all parameters; for convenience of comparison we also report here the results of the next sections. 

In Fig. Q instead we fix ft,, i.e. we assume that its fiducial value has a negligible error. If we fix h and we get 
Fig. ([5]); the other contours do not improve with respect to fixing h only (see Table W7\ . The significant improvement 
on the errors shows that there is much to gain by reducing the errors on ft,, Ug. The error on goes from 0.009 (fully 
marginalized) to 0.004 (fixing both h and n^). with a consequent better estimation of Vic- 
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Figure 2: CMB unlensed TT temperature spectra for three values of j3. Data are taken from WMAP7 (i^ . 



Parameter 


oi CMB 


P{k) 


a, WL 




0.0094 


0.0015 


0.012 


a 


0.55 


0.12 


0.083 


no 


0.022 


0.010 


0.012 


h 


0.15 


0.036 


0.039 




0.00087 


0.0022 


0.010 


Us 


0.014 


0.034 


0.026 


0-8 




0.0084 


0.024 


log^ 


0.0077 







Table TIL 1-cr errors for the set Q = a, ilc, /i, ^6, ?is, erg, log A} of cosmological parameters, using CMB data, P{k) and 
Weak Lensing (WL). Here and in all the following tables the errors are fully marginalized. 



Parameter Oj CMB a, CMB(fe) Oj CMB(fa,Ws) 





0.0094 


0.0044 


0.0041 


a 


0.55 


0.52 


0.44 




0.022 


0.011 


0.0079 


h 


0.15 








0.00087 


0.00081 


0.00043 


Us 


0.014 


0.012 




CT8 











log{A) 


0.0077 


0.0074 


0.0057 



Table IV: 1-a errors using CMB for the set © = a, flc, h, fit,, Us, as,\og{A)} (left column), Q = {/S'^ , a, flc, fit, ns, erg, log(yl)} 
(middle column, h fixed to the reference value), & = {P^ , a,Qc, h,Qb,o's,log{A)} (right column, h and Us both fixed to their 
reference values). 
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Figure 3: Predicted confidence contours for tlie cosmological parameter set Q = {/S^ , a, ^Ic, h,^lf,,ns, log A} using CMB Planck 
specifications. 

V. RESULTS FOR THE GALAXY POWER SPECTRUM 

The main effects of the coupling on the matter power spectrum are the shift in the matter-radiation equahty, the 
change in location and amplitude of the BAO and the speed up in the perturbation growth. As already noted, since 
dark matter dilutes faster for larger there is more dark matter in the past than in uncoupled cases and therefore 
the equality moves to higher redshifts. This implies that the wavelenghts at which perturbations reenter during 
radiation domination are smaller and thus the power spectrum turnaround moves to smaller scale. Just as for the 
CMB, the acoustic peaks also move to smaller scales and their amplitude is reduced. These effects are clearly visible 
in Figs. (|6l7p . The difference in the parameter growth function is plotted in Fig.®. In Fig.([7|) we make manifest the 
effect of the coupling on baryonic acoustic oscillations. For each /3 we plot the ratio between the power spectrum and 
its smooth spectrum (obtained reducing the amount of baryons in the CAMB code) . The ratios are then normalized 
to high momenta k. As we can see the coupling affects both the amplitude and position of the BAO peaks, shifting 
them to higher momenta. The effect of /? on P{k, z) is therefore quite strong since it occurs simultaneously on the 
wiggle position and amplitude, on the broad spectrum shape, and on the growth. This results in a tight constraint 
on the coupling parameter. 

Following j50i | , we build a function that provides the observed linear spectrum at all redshifts and for all cosmological 




Figure 4: Predicted confidence contours for tlie cosmological parameter set O = {0^ , a, Q,c, h, Q,h, Ws, log A] (in black, equal to 
Fig. ([3])) using CMB Planck specifications. Overplotted in dashed are the predicted confidence contours for the cosmological 
parameter set {/3^, a, Q,c, Ht, n^, log A}, with h fixed to the reference value. 



parameters: 



Pr.obs{kr: fJ-r'i — Ps{z) + 



D^{z)Hr{z) 



G\z)h\z)ai{l + Mz)fi'yP{k, z - 0) 



(18) 



where G{z) is the growth factor, b{z) is the bias, (3d{z) is the redshift-distortion factor, P{k, z = 0) is the undistorted 
Hnear matter spectrum at z = normalized to unity, fi is the direction cosine, D and H are the angular diameter 
distance and the Hubble rate at the shell redshift z, respectively, and Ps{z) is the z- dependent shot-noise correction, 
to be marginalized over in every redshift bin. The subscript r (for 'reference') indicates quantities calculated in the 
fiducial model. In linear theory we have (id{z) = f{z)/b{z) where f{z) — dlogG/dloga is the growth rate. 
In this case, the Fisher matrix for every redshift bin shell is (soj . 



+1 



d/x 



fc^dfc 



d In Fobs (fc, /i) 9 In Fobs {k,n) 



80, 



n{z)Fobsik,n) 



n{z)Fobs{k, /i) + 1 



(19) 
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Figure 5: Predicted confidence contours for the cosmological parameter set = {0^ ,a,Q,c,h,^b,ns,\ogA} (in black, equal 
to fig. (SI)) using CMB Planck specifications. Overplotted are the predicted confidence contours for the cosmological param- 
eter set {/3'^, a, f2c, rij,, ris, log A} (dashed contours, with h fixed to the reference values) and the cosmological parameter set 
a, f^c, fit, log A} (dotted contours, with h and Us fixed to the reference values). 

where n{z) is the galaxy number density at redshift z and Vs is the survey volume of the redshift shell. The power 
spectrum at 2; = and the functions G{z), H{z), D{z) are all obtained numerically by solving the background and 
perturbation equations of the system for each value of the parameters and the derivatives in Fij are evaluated 
numerically. We follow [5l| and instead of marginalizing over G and / we express them in function of the 

parameters; in particular we form the combination g(z) = Gfag for each redshift bin and write 



and we marginalize over Pd{z) assuming an independent parameter (id for every redshift bin. The prescription for 
kmax is that the variance in cells ci^ikmaxTz) = 0.25, resulting in a kmax ~ 0.16/i/Mpc for the first shell at z = 0.6. 
This conservative cut of the higher momenta allows us to discard the problem of the non-linear correction. 

Here and in the next section we have a new parameter, namely trgj whose fiducial is fixed to 0.8. In principle, the 
parameter log A employed in the CMB section is related to cts and therefore we should not count them separately. 
However, the two normalizations are taken at very different epochs and we are being conservative by assuming them 
to be independent. 

We use specifications for a Euclid-like mission based on the Euclid Red Book 3. The fiducial bias b{z) foUows f5^; 
we adopt redshift bins of Az — 0.1 between 0.5 < z < 2 and a sky coverage of 15,000 sq. deg. [3]. The expected 
number of Ha emitters per deg^ between z — dz/2 and z + dz/2 is computed for Ha fiux Fhu = 4.00 x IQ^^^erg/ s/cm? 
according to [5^. In Fig. ^ we show the expected confidence regions from galaxy power spectra. 

The Euclid- like P{k) probe gives in general better constraints than Planck (and also WL, see below), as pointed 
out in 0. This applies also to the coupling parameter /3: we find indeed < 0.0015, a constraint six times stronger 
than for Planck CMB. 



G''{z)h\z)al{\ + Pd{z)^^'' fP{k, z = 0) 




{l + Pd{z)^l^fP{k,z 



0) 



VI. RESULTS FOR WEAK LENSING 



An excellent probe for dark energy parameters that complements the aforementioned techniques is weak lensing (see 
[s^ l for an extensive treatment), due to the sensitivity with respect to the growth of structure and equal treatment of 
dark and baryonic matter. The weak lensing power spectrum depending on the multipole (. can be written in terms 
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Figure 6: The combination G{z)f{z) that appears as the amplitude of the power spectrum, for various values of /?. G{z) is 
normalized to unity today. 




Figure 7: BAO oscillations in P{k) for three values of /3. The curves show the ratio P{k)/ Ps{k), where Ps{k) is the corresponding 
smooth power spectrum. 



of the matter power spectrum as [5 



where 



is the window function and 



Wiz)^ I dz(l-^)r.(z) (21) 



n(z) = z2gxp(-(z/zo)3/2) (22) 
the normalized galaxy distribution function [56]. Note that zq is related to the median redshift via Zmcd ~ 1.412zo. 
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Figure 8: Predicted confidence contours for the cosmological parameter set = q, Sic, h, ^If,, Us, as} using the galaxy power 
spectrum. 



In order to extract as much information from the cosmic shear as possible, we will also employ the so-called weak 
lensing tomography (57| . where we form M = 5 redshift bins and compute the cross-correlation and power spectrum 
of the shear field. With this, eq. (^0)) becomes |57| 



9 r°° 



W,{z)W,{z)H\z)nl^{z) 



il + z) 



Prr 



7rr(z) 



(23) 



where Wi{z) now depends on ni{z), the galaxy distribution in the i-th redshift bin. The binned galaxy distribution 
is normalized to unity and has been convolved with a Gaussian to account for photometric redshift errors Az(l -I- z) , 
i.e. 



(24) 



i-th bin 



where z^ is the center of the i-th redshift bin and Ai a normalization factor. In Fig. ^ we show the convergence 
power spectrum for three values of (3. The behavior is not trivial since the spectra are a convolution of background 
and perturbation evolution. 
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Figure 9: The convergence power spectra Pn and P55 for three values of /3. 
The Fisher matrix for the weak lensing power spectrum is then given by (58| 

^^p-h^y 2^ 2 dOo, ^'^ 890 ™ 

with the covariance matrix 

where 7int = 0.22 in the shot noise term is the intrinsic galaxy ellipticity [s^ and 

n» = 3600 f— j ne/N 



(25) 



(26) 



(27) 



with ne being the total number of galaxies per arcmin^, assuming that the redshift bins have been chosen such that 
each contain the same amount of galaxies. Since we consider multipoles up to £max = 5000, we need to apply non- 
linear corrections to the matter power spectrum, for which we use the fitting formulae from Ref. (59| . It is clear that 
this is far from satisfactory since these non-linear corrections are calibrated through ACDM TV-body simulations and 
should not be used outside these cases. However at the moment there are no suitable analytical extensions to coupled 
models (for a recent effort in this direction see e.g. [4lj) so we cannot improve on this. Moreover, since our fiducial 
is indeed almost a ACDM cosmology we can assume that the non-linear correction does not introduce a large bias in 
our results. 

As before, we select a Euclid-like mission as our probe. We use the specifications listed in Table fVl as taken from 
the Euclid Red Book [s"]. In Fig. pO|) we show the marginalized confidence contours. The WL constraint on is 
0.012, weaker than either CMB or P{k). Interestingly, however, the constraint on the potential slope a is the strongest 
among the various probes. 



VII. COMBINING CMB, P(K) AND WL 



In this Section we finally combine all the results of the various probes. The combination of P{k) and WL is 
not totally accurate since the two probes are assumed to be independent while in fact they observe the same field 
of galaxies and the lensing signal is correlated with the density distribution. However, since as we have seen the 
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Parameter 


Value 


Description 


A 


15 000 deg2 


Survey area 


ng 


30 


Galaxies per arcmin^ 


f 

^max 


5000 


Maximum multipole 


^mcd 


0.9 


Median redshift 


A, 


0.05 


Photometric redshift error 



Table V: Mission goals for a Euclid-like project. Q 




Figure 10: Predicted confidence contours for the cosmological parameter set Q = a, f2c, /i, ?is, crs} using weak lensing 
Euclid-like specifications. 



constraints from P{k) are in general quite stronger, we expect that an improved treatment will not change drastically 
our results. 

The combined Fisher confidence regions are plotted in Fig. pTj) and ([T^ and the results are in Table (|VII[) . The 
main result is that future surveys can constrain the coupling of dark energy to dark matter to less than 3 • 10~^. 

We can ask whether a better knowledge of the parameters {a, i^c, /i, ^^fj, f^-s, ctsi log(^)}, obtained by independent 
future observations, can give us better constraints on the coupling In Table IVlIl we list the l-cr error on obtained 
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Parameter 


oi CMB+P(ifc) 


Or CMB+P(A:)+WL 




0.00051 


0.00032 


a 


0.055 


0.032 


fie 


0.0037 


0.0010 


h 


0.0080 


0.0048 


fib 


0.00047 


0.00041 


Us 


0.0057 


0.0049 




0.0049 


0.0036 


\og(A) 


0.0051 


0.0027 



Table VI: l-cr errors for the set O = {/3^, a, Sic, /i, SI;,, rij, crs, log(yl)} of cosmological parameters, combining CMB+P(fc) (left 
column) and CMB+P(A;)+WL (right column). 




-2-1012 _2 -1 1 2 



Figure 11: Comparison among predicted confidence contours for the cosmological parameter set G = 
, Tis, (78, /o5r(yl)} using CMB (Planck, blue contours), P(A:) (pink-violet contours) and weak lensing (orange-red 
contours) with Euclid-like specifications. 
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0.30 



0.15 




0.0004 0.0008 0.0012 



a 0.226 




0.0004 0.0008 0.0012 




0.0004 0.0008 0.0012 





a 0.226 



0.224 



0.0004 0.0008 0.0012 o.222 




0.10 0.15 0.20 0.25 0.30 



0.0465 
0.0460 
0.0455 
'^0.0450 
0.0445 
0.0440 





0.10 0.15 0.20 0.25 0.30 



0.10 0.15 0.20 0.25 0.30 



0.720 
0.715 
0.710 
0.705 
0.700 
0.695 
0.690 
0.685 




Figure 12: Combined predicted confidence contours for the cosmological parameter set Q = {13^ , a,flc, h,flf,,ns, as,log{A)} 
from CMB (Planck), P{k) and weak lensing (Euclid- like) specifications. 



when we fix, one after the other, each of the remaining parameters to their reference value instead of marginaUzing. Pi 
corresponds to fixing one parameter (a) and marginalize over all the others, P2 corresponds to fixing two parameters 
(a and flc) and marginalize over all the others and so on until Pj where we estimate the error on when we fix all 
parameters. 

The first column corresponds to the effect on CMB only. CMB alone gains the most from a better knowledge of 
the parameters. In this case the error on improves by almost one order of magnitude when flc is fixed (thanks 
to the degeneracy with this parameter) and of almost two orders of magnitude when also h and fif, are known. A 
better knowledge of a, Sic, /i, Jlfc makes CMB constraints alone comparable to the combination of all probes (when all 
parameters are marginalized over). 

P{k) observations (second column) improve bounds by almost an order of magnitude when a, Sic, h, Qi, have been 
fixed, becoming even better than the ones obtained in the marginalized combination of all probes. The third column 
shows weak lensing constraints, which gain one order of magnitude when all parameters are known: note that in this 
case the amplitude of perturbations erg is still relevant even when all previous parameters have already been fixed. 

The last column combines all three probes. Interestingly, the combination of CMB, power spectrum and weak 
lensing is already very good when all parameters are marginalized over; no much gain is obtained, overall, by a better 
knowledge of the other parameters. This is confirmed also when looking at Table IVnll Here we show the errors on 
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Set of fixed parameters 


CMB 


P(fc) 


WL 


CMB + P{k) - 


F WL 


(Marginalized on all params) 


0.0094 


0.0015 


0.012 


0.00032 




A = {a} 


0.0093 


0.00085 


0.0098 


0.00030 




P2 = {Pi} U {Q4 


0.0026 


0.00065 


0.0084 


0.00030 




Ps = {P2} U {h} 


0.00054 


0.00040 


0.0076 


0.00026 




P4 = {PalUl^b} 


0.00033 


0.00028 


0.0037 


0.00020 




P5 = {P4}U K} 


0.00033 


0.00024 


0.0034 


0.00019 




Pe = {Ps} U {erg} 


0.00033 


0.00024 


0.0017 


0.00019 




P7 = {Ps} U {log(A)} 


0.00032 


0.00024 


0.0017 


0.00019 




Table VII: 1-a errors for /3^ for CMB, P(fc), WL and CMB+P(fc)+WL. 


For each line the set of parameters Pi has been fixed 


to the reference value. The first line corresponds 


to the case in which we 


have marginalized 


over all parameters. 




Fixed parameter 


CMB 


P{k) 


WL 


CMB + P{k) - 


F WL 


(Marginalized on all params) 


0.0094 


0.0015 


0.012 


0.00032 




a 


0.0093 


0.00085 


0.0098 


0.00030 






0.0026 


0.00066 


0.0093 


0.00032 




h 


0.0044 


0.0013 


0.011 


0.00032 




nt 


0.0087 


0.0014 


0.012 


0.00030 




Us 


0.0074 


0.0014 


0.012 


0.00028 






0.0094 


0.00084 


0.0053 


0.00030 




log(A) 


0.0090 


0.0015 


0.012 


0.00032 





Table VIII: 1-a errors for for CMB, P(fc), WL and CMB+P(fc)+WL. For each line, only the parameter in the left column 
has been fixed to the reference value. The first line corresponds to the case in which we have marginalized over all parameters. 



when we have a better knowledge of only one other parameter, which is here fixed to the reference value. All 
remaining parameters are marginalized over. The combination of CMB, power spectrum and weak lensing is already 
a powerful tool and a better knowledge of one parameter doesn't improve much the constraints on /3^. CMB alone, 
instead, improves by a factor 3 when flc is known and by a factor 2 when h is known. The power spectrum is mostly 
influenced by fie, which allows to improve constraints on the coupling by more than a factor 2. Weak lensing gains 
the most by a better knowledge of erg . 

The constraint on can be easily converted into a constraint on the Brans-Dicke coupling parameter ujbd since 
(see e.g. ^) 

3 + = ^ (28) 

We obtain then 

UJBD > 800 (29) 

The post-Newtonian parameter jppN is related (in the limit 7 ^ 1) to lubd as ^ppN = 1 — 2/(2a;B£i + 3) so that our 
final predicted constraint on 0^ reflects into a constraint 

I7PPJV - 1| < 1.2 X 10"^ (30) 

This is still two orders of magnitude weaker than the limits on a scalar coup lin g to baryons obtained in laboratory 
or solar system experiments, of the order of \^ppN — 1| < 10~^ (see e.g. [24l. [25l.l60j ) but of course the cosmological 
bounds are complementary to the local measurements since here we deal with the coupling to dark matter. 

VIII. CONCLUSIONS 

Future cosmological surveys will measure the properties of our universe to a unprecedented precision. This will 
allow not only to measure the main cosmological parameters but also to test models that challenge the standard 
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scenario of particle physics and cosmology. The two main tools of observational cosmology, namely cosmic microwave 
background and large-scale structure, can be combined to extract the maximal amount of information. 

In this paper we combined forecasts from the Planck CMB experiment and from future redshift surveys and weak 
lensing surveys based on the proposed Euclid satellite to get constraint on a model of coupled dark energy. In this 
model, dark energy mediates a force acting on dark matter and produces an additional attractive interaction which 
modifies Einstein's gravity. The coupling parameter measures the amount of this modification. We find that 
combining Planck with a Euclid-like survey can constrain to a level of 3 • 10~^, two orders of magnitude better than 
current constraints [20, 22]; this bound is weaker but complementary to the small-scale limits set by local gravity tests 
on Yukawa corrections (see e.g. (24l. |25|V A better knowledge of flc, h, Ug can improve the precision of measurement 
of the coupling as obtained by CMB measurements, due to manifest degeneracies among these parameters and the 
fifth force as measured by 0^ . Interestingly though, the combination of CMB, power spectrum and weak lensing is 
already so powerful that no much gain is obtained, overall, by a better knowledge of the other parameters. In this 
sense, these three probes are optimally complementary probes of the dark energy-dark matter coupling. 
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